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Abstract — In this article, we develop a continuous implicit seventh-eight method using interpolation and 
collocation of the approximate solution for the solution ofy’ = f(x, y) with a constant step-size. The method 
uses power series as the approximate solution in the derivation of the method. The independent solution was 
then derived by adopting block integrator. The properties of the method was investigated and found to be 
zero stable, consistent and convergent. The integrator was tested on numerical examples ranging from linear 
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solution is closer to the exact solution and also, the absolutes errors perform better than the existing methods. 
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I. INTRODUCTION 

Numerical analysis predates the invention of modern 
computers by many centuries. Linear interpolation was 
already in use more than 2000 years ago. Many great 
mathematicians of the past were preoccupied by numerical 
analysis, and this is obvious from the names of important 
algorithms like Newton's method, Lagrange interpolation 
polynomial, Gaussian elimination and Euler's method. 
Despite the long and rich history of numerical analysis, 
"modern" numerical analysis is characterized by the 
synergy of the programmable electronic computer, 
mathematical analysis, and the opportunity and need to 
solve large and complex problems in applications. The need 
to advance in applications such as ballistics prediction, 
neutron transport, and non-steady multidimensional fluid 
dynamics motivated the development of the computer 
which depended strongly on advances in numerical analysis 
and mathematical modelling. 

1.1 Aim and Objectives 


The aim of this research work is to develop continuous 
method via fractional-step methods by using power series 
polynomial as basic functions for solving IVPs in ODEs. 
The objectives are to. 

i. develop hybrid block methods through continuous 
implicit seventh-eight method 

ii. analyzing the basic properties of the method; zero 
stability, order, convergence and region of absolute 
stability 

iii. implementing a self-staring implicit methods 

without the rigor of developing predictors 

iv. developing computer codes to implement the 

method derived using MATLAB 

v. examine the exact solution and the computed 

solution, if they are in agreement 

vi. compute the time taken for computation in seconds 

1.2 Significance of the Study 

The study is significant for the following reasons 

1. This work centered on the numerical solution of 
first order ordinary differential equations of initial 
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value problem. Hence, the differential equations 
whose order is more than one are not considered. 

2. The proposed method makes use of power series 
polynomial, other polynomials are not considered 
in this work. 

3. Throughput this research work, linear fractional- 
step method is considered, Non hybrid points are 
not considered. 

4. Most of the modeled problems in ordinary 
differential equations do not have analytical or 
theoretical solution, hence the use of numerical 
solution is very important 

n. LITERATURE REVIEW 

This work is based on numerical solution of first order 
initial value problems of ordinary differential equation, 
which is of the form 

y'= f(x,y), y(a) = T], a < x < b (1) 

where f is continuous with the interval of integration. It is 
assumed that f satisfies Litchitz condition which guarantees 
the existence and uniqueness of solution (1). 

Countless real life problems in sciences, engineering, 
biology and social sciences are model of first order ordinary 
differential equations. For instance, Environmental Models 
such as SIR Model, Growth model, Mixture model just to 
mention few. Olver (2008) and Kandasamy et al. (2005), 
reported that equation (1) is used in simulating the growth 
of populations, trajectory of a particle, simple harmonic 
motion, deflection of a beam etc. With the advent of the 
modern high speed electronic digital computers, the 
numerical integrators have been successfully applied to 
study problems in mathematics, engineering, computer 
science and physical sciences such as physics, biophysics, 
atmospheric sciences and geosciences, Jain et al. (2007) 
Many numerical integration schemes to generate the 
numerical solution to problems of the form (1) have been 
proposed by several authors. The one-step methods include 
those developed by Fatunla (1994), Aashikpelokhai (1991), 
Ayinde and Ibijola (2015), and Kama and Ibijola (2000), 
Ibijola (1998), Ogunrinde (2010) etc . These methods were 
constructed by representing the theoretical solution y(x) to 
(1) in the interval, [ X n , X n+l ] Tl > 0 by linear and non-linear 
polynomial interpolating functions. 


Furthermore, authors like Butcher (2003), Sunday et al. 
(2015), Lukman and Olaoluwa, Areo et al. (2011), Ibijola, 
Skwame and Kumleng (2011), Olangegan et al, (2015), 
Areo and Omojola (2017), Adeniran and Longe (2019), 
Adeniyi et al. (2006), Onumanyi et al, (1994), Sunday and 
Odekunle (2012), Badmus (2014), Jator (2014), Kayode and 
Adeyeye (2013) and Ibijola et al. (201) have all proposed 
linear multistep methods (LMMs) to generate numerical 
solution to (1). These authors proposed methods in which 
the approximate solution ranges from power series, 
Chebychev’s, Lagrange’s, Lucas polynomial, Legendre 
Polynomial and Laguerre’s polynomials. In spite of the 
improvement made by diverse researchers in developing 
schemes that solve (1), the order of the schemes are low 
resulting in lower accuracy when used to solve problems 
which is a setback. Hence, the setbacks by previous works 
motivated a need to develop a seventh-eight method for 
solving first order ordinary differential equations unlike 
other methods which is one-step or multi-step. The method 
will be analysed and test on some numerical examples to 
test the accuracy and efficiency of the new method in order 
to ascertain how reliable it is. 


III. MATERIALS AND METHODS 

The power series of the form (2) is considered 

5+a-l 

y(x)= Yj a „Pn ( x )> 

n=0 

( 2 ) 

The first derivative of (2) is given by (3) 


5+o -1 

y'(.X)= Y a nP'n( X )’ 

(3) 

as an approximate solution to equation (1) where 
p n (x) = x n ,. 8 and a are the number of distinct 
collocation and interpolation points respectively. 
Substituting the first derivative of (3) into (1) gives 


f(X,y(X%)= Y a nP'n ( X ). 


( 4 ) 


http://www.aipublications.com/ijcmp/ 





Open Access 


International journal of Chemistry, Mathematics and Physics (IJCMP) 
[Vol-4, Issue-3, May-Jun, 2020] 
h ttps://dx. doi. org/10.22161/ijcmp. 4.3.2 

ISSN: 2456-866X 


nl 1 ) 7 3 

Collocating (4) at X — U — — and interpolating (3) at X — — leads to a system of nine equations which is solved with the aid 

V 8 7 8 4 

of Maple Mathematical Software version 18.0 to obtain a n , n = 0, 1, • • •, 8 . The Cl n ' s obtained are then substituted into (3) to 
obtain the continuous form of the method 

mf n +A(0/ +i + A (0/ +i + A (0/ +3 + A (0/ +i +A (0/ + 

+A (Of 3+A (f)f v 


y(x) = a 3 (t)y 3 +h 


where OC n and J3 n are continuous coefficients. 

a (0 = 1 


P 0 (0 = ~ 3Q ^ 4Q (6144 A 4 t 4 - 12288 A 3 t 3 + 8192 A 2 t 2 - 2048 A t 
+ 123) (4ht — 1 ) 2 (4 A / — 3 ) 2 

jSj (0=3^- (86016 A 6 / 6 -202752 A 5 / 5 + 176128 A 4 / 4 - 69312 A 3 f* 

y 

+ 11328 A 2 ^ 2 — 216 A ? — 81) (4A/-3) 2 

P. (0 = - (229376 A 6 / 6 - 507904 A 5 / 5 + 399360 A 4 f 4 - 133120 A 3 f* 

_L 3360 

4 

+ 15728 A 2 ; 2 + 24A/ + 9) (4A/-3) 2 
P 3 (t)= -j^- (215040 A 6 / 6 - 445440 A 5 f 5 + 317440 A 4 / 4 - 92448 A 3 f* 

T 

+ 9528 A 2 t 2 — 136 A / — 51) (4A/-3) 2 

p x (0 = - 3Q j, 4Q (3440640 A 6 / 6 - 6635520 A 5 t 5 + 4311040 A 4 / 4 
2~ 

- 1155072 A 3 / 3 + 118032 A 2 / 2 + 216 A / + 81) (4A/-3) 2 
j3 5 (0 = (28672 A 6 1 6 - 51200 A 5 / 5 + 30720 A 4 / 4 - 8000 A 3 f 5 

y 

+ 736 h 2 1 2 —24ht — 9) (4A/-3) 2 
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J M ?) = 'l0240' (4ht-3) (2752512 h 1 1* — 6586368 h 6 1 6 

+ 5955584 h 5 t 5 - 2586624 h 4 1 4 + 545088 h 3 t 3 - 47696 h 2 1 2 - 492 ht 

- 369) 

P 7 (t)= t 2 h 2 (768 h 4 t 4 - 1152 h 3 t 3 + 640 h 2 t 2 - 156 h t + 15) (4 h t 

¥ 

— 3) 2 

The continuous method (5) is used to generate the main method. That is, we evaluate at X = X 7 

y „ + L =y „ + L + ~9^M h { 1375/ » + 123133/ „ + .L + 41499/ „ + L 

8 4 { + 2 + 4 

-11351/ j + 139849/ 3 - 88547/ 3 - 121797/ 5 

« + — n + — n + — n + — 

+ 36799/ +1 j 


Continuous method (5) is also used to generate the additional method at the non-interpolation points. That is, we evaluate at 
X = X J,X £,X J and JC = X 5 


T„=T 3_~Ym h ( 4l f n + 21 f _ L + 27 / j_+216/ J _+41/ ^ 

n+ 4 { H+ 2 H+ 4 H+ 8 n+ 4 

+ 272/ 3_+216/ 

" + 8 ” + 8 J 

y „ + ±= y „ + ± + id * h f 55/ » - W5/ . + ± - 5805/ „ + J_ 

8 4 2 4 

-1871/ ^ - 1871/ ^ - 4475/ ^ - 5805/ ± + 55/ ^ 

«+ g ” + 4 «+ g " + 8 " + gj 

^ + ±+ + i-lW/ 33 / + ± + 17 / + ±-/ + ± + 17 / + i 

4 4 ( 2 4 8 + 4 

+ 612/ 3_+612/ ^-4 / ^ 

” + 8 ” + 8 " + 8 J 


(7) 


( 8 ) 


( 9 ) 
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8 4 1, 2 4 

- 1685 / 3 _- 2777 / ^ - 5535 / ^ + 45 / ^ 

n + 4 H+ g «+ g W+ 8 J 
^ + ±^ + A + ^k/v»- 115 / + i + 13 / + ±- 4 /_ 

2^4 l, 2 4 

- 1363 / 3 - 208 / 3 - 4968 / 5 + 32 / ? \ 

n + J n + J n + J n + J) 

3 ' B+ i^„ + i + %4o/ 351/ » + 44797/ „ + i + U547/ „ + ± 

8 4 2 4 


- 2999/ , 

” + 8 

+ 1375/ ? 


t 2 ” 4 

-47799/ 3 -26883/ 3 - 101349/ 5 

« + — « + — ” + T 


The methods derived in equation (6) to (12) are combined and implemented as a block in solving numerical problems of (1). 


IY. ANALYSIS OF THE METHOD 


4.1 Order and Error Constant of the Block Method 


We adopt the method proposed by Fatunla (1994), Lambert (1973) and Henrici (1962) to obtain the order of the method obtained 
from (6-12) as (8,8,8,8,8,8,8) T and Error constants as 


[-6.97e“ u , -5.47e“ u , -6. 14e~ n , -5.62e~ n , -6.29e n , -4.79<T n , -1. 1 Se~ n J 


(13) 


4.2 Zero Stability of the Block Method 

A block method is said to be zero stable, if the roots 


det [AA (0) - A® ] = 0 


(14) 


of the first characteristic polynomial satisfy |/l| < 1 and for the roots with | < 1, the multiplicity must not exceed the order of 

the differential equation. 


http://www.aipublications.com/ijcmp/ 









Open Access 


International journal of Chemistry, Mathematics and Physics (IJCMP) 
[Vol-4, Issue-3, May-Jun, 2020] 
h ttps://dx. doi. org/10.22161/ijcmp. 4.3.2 

ISSN: 2456-866X 



'1 

0 

0 

0 

0 

0 

0" 


'0 

0 

0 

0 

0 

o r 



0 

1 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 1 



0 

0 

1 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 1 


z 

0 

0 

0 

1 

0 

0 

0 

- 

0 

0 

0 

0 

0 

0 1 



0 

0 

0 

0 

1 

0 

0 


0 

0 

0 

0 

0 

0 1 



0 

0 

0 

0 

0 

1 

0 


0 

0 

0 

0 

0 

0 1 



v 0 

0 

0 

0 

0 

0 

K 


v0 

0 

0 

0 

0 

0 1, 



A = Z 6 (z -1) = 0, Z = o, 0,0,0,0,0,0 
Hence the block is zero stable. 

4.3 Consistency of the Method 

Our new developed block method is consistent since the order of each of the method is greater than 1. 

4.4 Convergence 

According to Lambert, 1973, Brown 1977 and Chollom et al. (2007) the necessary and sufficient condition for a linear multistep 
method to be convergent is for it to be consistent and zero stable. 

Hence Convergence = Consistency + Zero-stability. 

Hence the new block method is convergent 

4.5 Region of Absolute Stability of the Proposed Method 

According to Ibijola et al. (2011), Areo and Omojola (2017), the stability matrix is expressed as 


M(z) = V + zB(M-zAy l U 

together with the stability function 


p(rj, z) = det( 7 1-M(z)) 

We express the block method obtained from (6-12) in the form 



(16) 


(17) 


(18) 


where 
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(19) 


The elements of the matrices A, B, U and V are substituted and computing the stability function with Maple software yield the 
stability polynomial of the method which is then plotted in MATLAB environment with the Newton-Raphson Method where 
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N — 400 and tol= 10 12 to produce the required absolute stability region of the method as shown in the figure 1 below. The 
graph is shows the region in which the method is absolutely stable. We conclude that the method is P-stable since the interval of 
periodicity lies between (0, 2.0). 


Region of Absolute Stability of the Block Method 



Y. NUMERICAL EXPERIMENT 
This section deals with the implementation of the block 
method in solving initial value problems (IVP) of first order 
ordinary differential equations. The method is coded in 
MATLAB (R2012a) version environment using window 10 
as an operating system on Acer Laptop. The new developed 
method is tested on some problems to determine the 
performance of the new proposed methods. The error is 
defined as Error = |j(x) — )> n (x)|, where y(x) is the 
exact solution, y n (x) is the computed result. 

5.1 Numerical Examples 

The following problems are chosen as test problems 

Problem 1: Consider the first order problem 


y' m ^y t y(G)= 1, 0 < X < 1, * = 0.1, 

( 20 ) 

Exact solution: y(x) = e ' 

Problem 2: Consider the Prothero-Robinson Oscillatory 
ODE solved by Sunday et al. (2015) 

y ' = L(y - sin x) + cos x, y(0) = 0, L = -1, * = 0.1 
( 21 ) 

Exact solution: y(x) = sin X 
Problem 3: Consider the Growth Model 
A bacteria culture is known to grow at a rate proportional to 
the amount present. After one 

hour, 1000 strands of the bacteria are observed in the 
culture; and after four hours, 3000 strands. 
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Find the number of strands of the bacteria present in the 
culture at time t : 0 < t < 1 


Let N(t) denote the number of bacteria strands in the culture 

at time t, the initial value 

problem modeling this problem is given by, 


-= 0.366jV, N(0) = 694 

dt 


dS[ 
dt 
<U_ 
dt 
dR 
dt 


= ju(l-S)-/3IS 
= -jil -yl + piS 
= -/jR + yl 


(23) 


( 22 ) 

The exact solution is given by 

N{t) = 694e 0366r 

Problem 4: Consider the SIR Model 

The SIR model is an epidemiological model that computes 

the theoretical number of people infected with a contagious 

illness in a closed population over time. 

The name of this class of models derives from the fact that 
they involve coupled equations relating the number of 
susceptible people S(t), number of people infected I(t) and 
the number of people who have recovered R(t). This is a 
good and simple model for many infectious diseases 
Including measles, mumps and rubella. It is given by the 
following three coupled equations: 


where fl,y,(f are positive parameters to be determined. 

Define y to be, 

y=S+I+R 

and adding the equations in (23) above, we obtain the 
following evolution equation for y, 

y' = jui\-y) 

Taking /J = 0.5 and attaching an initial condition 
y(0) = 0.5 (for a particular closed population), we obtain, 

-j- = 0.5(1 — y), y(0) = 0.5 

dt 

whose Exact solution is 


y(t) = 1 —0.5e“° 5f 


Table 1: The exact solution, computed solution, the absolute error and the time taken for the computation in (seconds) of the new 
method for problem 1 


X-value Exact-solution Computed-solution Error Time/sec 


0.10 

0.20 

0.30 

0.40 

0.50 

0.60 

0.70 

0.80 

0.90 

1.00 


0.9048374180359595 

0.8187307530779814 

0.7408182206817172 

0.6703200460356384 

0.6065306597126323 

0.5488116360940253 

0.4965853037914083 

0.4493289641172203 

0.4065696597405978 

0.3678794411714410 


0.9048374164190154 

0.8187307469199735 

0.7408182076888079 

0.6703200243242226 

0.6065306278521965 

0.5488115930362602 

0.4965852486351736 

0.4493288965201419 

0.4065695792037555 

0.3678793476848027 


1.6169441e-09 
6.1580079e-09 
1.2992909e-08 
2.1711416e-08 
3.1860436e-08 
4.3057765e-08 
5.5156235e-08 
6.7597078e-08 
8.0536842e-08 
9.3486638e-08 


0.205696 

0.367746 

0.551484 

0.702339 

0.872062 

1.089429 

1.311644 

1.510724 

1.704233 

1.928272 
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Table 2: The exact solution, computed solution, the absolute error and the time taken for the computation in (seconds) of the new 
method for problem 2 


X-value Exact-solution Computed-solution Error Time/sec 


0.10 

0.20 

0.30 

0.40 

0.50 

0.60 

0.70 

0.80 

0.90 

1.00 


0.0998334166468282 

0.1986693307950614 

0.2955202066613387 

0.3894183423086477 

0.4794255386041985 

0.5646424733950294 

0.6442176872376839 

0.7173560908995147 

0.7833269096274749 

0.8414709848078880 


0.0998334154028488 

0.1986693259603136 

0.2955201961494993 

0.3894183242933311 

0.4794255115178666 

0.5646424359271566 

0.6442176383320178 

0.7173560297503054 

0.7833268356745610 

0.8414708977305652 


1.2439794e-09 
4.8347478e-09 
1.0511839e-08 
1.8015317e-08 
2.7086332e-08 
3.7467873e-08 
4.8905666e-08 
6.1149209e-08 
7.3952914e-08 
8.7077323e-08 


0.181433 

0.380374 

0.564763 

0.751144 

0.952488 

1.129044 

1.332819 

1.620844 

1.904968 

2.125000 


Table 3: The exact solution, computed solution, the absolute error and the time taken for the computation in (seconds) of the new 
method for problem 3 


X-value Exact-solution Computed-solution Error Time/sec 


0.010 

0.020 

0.030 

0.040 

0.050 

0.060 

0.070 

0.080 

0.090 

0.100 


696.5446939492859400 

699.0987185430900600 

701.6621079941421600 

704.2348966406204900 

706.8171189466110000 

709.4088095025695100 

712.0100030257847300 

714.6207343608434700 

717.2410384800974700 

719.8709504841319800 


696.5446939492180700 

699.0987185427908300 

701.6621079934537900 

704.2348966393841500 

706.8171189446453600 

709.4088094997416600 

712.0100030218745800 

714.6207343557549600 

717.2410384736124300 

719.8709504760998900 


6.7871042e-011 0.022520 
2.9922376e-010 0.050186 
6.8837380e-010 0.070106 
1.2363444e-009 0.090180 
1.9656454e-009 0.110484 
2.8278464e-009 0.133450 
3.9101451e-009 0.152786 
5.0885092e-009 0.175301 
6.4850383e-009 0.210087 
8.0320888e-009 0.232457 


Table 4: The exact solution, computed solution, the absolute error and the time taken for the computation in (seconds) of the new 
method for problem 4 


X-value Exact-solution Computed-solution Error Time/sec 


0.010 

0.020 

0.030 

0.040 

0.050 


0.5024937604036588 0.5024937604048754 1.2165824e-012 0.043527 
0.5049750831254160 0.5049750831324521 7.0361494e-012 0.048093 
0.5074440301984686 0.5074440302153604 1.6891821e-011 0.053913 
0.5099006633466223 0.5099006633774158 3.0793479e-011 0.059570 
0.5123450439858337 0.5123450440363059 5.0472182e-011 0.063933 
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0.060 0.5147772332257459 0.5147772332973700 7.1624151e-011 0.080116 

0.070 0.5171972918712168 0.5171972919729365 1.0171974e-010 0.085281 

0.080 0.5196052804238385 0.5196052805535286 1.2969015e-010 0.093241 

0.090 0.5220012590834500 0.5220012592496057 1.6615576e-010 0.097912 

0.100 0.5243852877496430 0.5243852879546123 2.0496926e-010 0.104638 


VI. SUMMARY AND CONCLUSION 
This paper examined the derivation, analysis and has 
demonstrated a successful implementation of seventh-eight 
step in form of linear multistep method. The method uses 
block approach for the solution of first order differential 

/ x . _ 1 1 3 1 5 

equations. The results for y I X n+j 1,7 = U, — 

8 4 8 2 8 

and — were obtained in block forms. The new method 

8 

speed up the computational processes, less burden in the 
implementations and also increase the rate of convergence 
of the solutions. The running times for the problem solved 
are less than four seconds as shown in Table 1-4. The exact 
solution of the method is in good agreement with the 
numerical solution. The method was found to be zero stable, 
consistence, convergent and P-stable as show in figure 1. 
Hence we conclude that seventh-eight method is reliable 
and gives good results with lesser time for computation of 
first order ordinary differential equations. 

REFERENCES 

[1] Olver P. J. (2008): Lecture Notes on Numerical Analysis, 
May, 2008, 

http://www. users.math.umn.edu/--olver/num.html 

[2] P. Kandasamy, K. Thilagavathy and K. Gunavathy, 
Numerical Methods, S. Chand and Company, New Delhi- 
110055, 2005. 

[3] Jain M. K, Iyengar S. R. K, and Jain R. K. “Numerical 
methods for scientific and engineering computation”. 5th 
Edition, New Age International, New Delhi, India, RK 
(2007). 

[4] Fatunla, S. O, Block Method for Second Order Initial Value 
Problem. International Journal of Computer Mathematics, 
England. Vol. 4, 1991,pp 55-63. 

[5] U.S.U. Aashikpelokhai, A class of nonlinear one-step rational 
integrator, PhD Thesis, University of Benin, Nigeria, 1991. 


[6] Ibijola, E. A, Skwame, Y. and Kumleng, G. (2011). 
Formation of hybrid block method of higher step sizes, 
through the continuous multi-step collocation, Am. Journal 
Sci. Ind. Res., 2(2), 161-173 

[7] Onimanyi, P., Awoyemi, D. O., Jator, S. N. & Sirisena, U. W. 
(1994). New Linear Multistep Methods with continuous 
coefficients for first order initial value problems. J. Nig. 
Math. Soc. 

[8] Adeniyi, R. B. Adeyefa, E. O. and Alabi, M.O. (2006), A 
continuous formulation of some classical initial value solvers 
by non-perturbed multistep collocation approach using 
chebyshev Polynomial as basis function. Journal of the 
Nigerian Association of Mathematical Physics, 10:261-274. 

[9] Brown, R. L. (1977): Some characteristics of implicit 
multistep Multi derivative Integration formulas. SIAM, J. on 
Numerical Analysis. 14, 982 — 993 

[10] Henrici, P. (1962), Discrete variable methods for ODEs. John 
Wiley, New York. Sunday J. and M.R. Odekunle (2012): A 
New Numerical Integrator for the Solution of Initial Value 
Problems in Ordinary Differential Equations. Pacific Journal 
of Science and Technology, Vol. 13 NO. 1, 2012, pp. 221- 
227. 

[11] Lambert, J. D. (1973): Computational methods in ordinary 
differential equation, John Wiley & Sons Inc. 

[12] Chollom, J. P„ Ndarn, J. N. and Kumleng, G. M. (2007): On 
some properties of the block linear multistep methods. 
Science World Journal, 2(3), 11 - 17. 

[13] Butcher, J. C. (2003). Numerical Methods for Ordinary 
differential systems. West Sussex, England: John Wiley & 

[14] Shampine, L.F. and Watts, H.A. (1969), Block implicit one- 
step methods. Journal of Computer Maths. 23:731-740. 

[15] Fatunla, S.O. (1994), A class of block methods for second 
order IVPs, Int. J. Comput. Math., 55: 119-133. 

[16] Areo E. A and Omojola M. T. (2015): A new One-Twelfth 
Continous Block method for the solution of modelled 
problems of ordinary differential equations, American Journal 
of Computational Mathematics, 5, 447-457 

[17] Olanegan, O. 0.,Ogunware, B. G., Omole E. O. Oyinloye, T. 
S. and Enoch,B. T.(2015): some Variable Hybrids Linear 
Multistep Methods for Solving First Order Ordinary 

[18] Differential Equations Using Taylor’s Series, IOSR Journal 
of Mathematics, 5 (1), PP 08-13 


http://www.aipublications.com/ijcmp/ 






Open Access 


International journal of Chemistry, Mathematics and Physics (IJCMP) 
[Vol-4, Issue-3, May-Jun, 2020] 
h ttps://dx. doi. org/10.22161/ijcmp. 4.3.2 

ISSN: 2456-866X 

[19] Adeniran A. O. & Longe I. O. “Solving directly second order 
initial value problems with Lucas polynomial. Journal of 
Advances in Mathematics and Computer science, 32(4), 1-7, 

(2019). http://dx.doi.org/10.9734/iamcs/2019/v32i430152 . 

[20] Lukman M. A, Olaoluwa O. E. “A One-step extended block 
hybrid formulae for solving Orbital problems”, Moj App. 

Bionics & Biomechanics,3(l):37-43, (2019). 

http://dx.doi.org/10.15406/moiabb.2019.03.0Q101 . 

[21] Sunday J., Skwame Y., and Huoma I. U., “Implicit One-Step 
Legendre Polynomial Hybrid Block Method for the Solution 
of First-Order Stiff Differential Equations”. B. J.of Maths. & 

Comp. Sci. 8(6), 482-491, (2015). 

https://doi.org/10.9734/BJMCS/2015/16252 

[22] Badmus A. M. “A New Eighth Order Implicit Block 

Algorithms for the Direct Solution of Second Order Ordinary 
Differential Equations”. American Journal of Computational 
Mathematics, 4, 376-386, (2014). 

http://dx.doi.org/10.4236/aicm.2014.44032 

[23] Jator, S.N., Lee, L. “Implementing a seventh-order linear 
multistep method in a predictor-corrector mode or block 
mode: which is more efficient for the general second order 
initial value problem”. SpringerPlus 3(1), 1-8 (2014). 
https://doi.org/10.1186/2193-1801-3-447 

[24] Adeyeye, O., Kayode, S.J. “ Two-step two-point hybrid 
methods for general second order differential equations” . 

Afr. J. Math. Comput. Sci. Res. 6(10), 191-196 (2013). 
https://doi.org/10.5897/AJMCSR2013.0502 

[25] Ayinde S. O., Ibijola E. A.. A New Numerical Method for 
Solving First Order Differential Equations. American Journal 
of Applied Mathematics and Statistics. 2015; 3(4): 156-160. 
doi: 10.12691/ajams-3-4-4 

[26] Kama, P. and Ibijola, E. A. (2000). On a New One - Step 
Method for Numerical Integration of Ordinary Differential 
Equtions. International Journal of Computer Mathematics, 

Vol.78, No. 3,4. 

[27] E.A.Ibijola, On the convergence consistency and stability of a 
one step method for integration of ordinary differential 
equation. International Journal of Computer Mathematics, 73 
(1998) 261-277. 

[28] R.B.Ogunrinde, A new numerical scheme for the solution of 
initial value problem (IVP) in ordinary differential equations, 

Ph.D. Thesis, Ekiti State University, Ado Ekiti. (2010) 


http://www.aipublications.com/iicmp/ 











